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I. INTRODUCTION 


Models of opinion formation have been studied by physicists since the 80’s and are now 
part of the new branch of physics called sociophysics. This recent research area uses tools and 
concepts of statistical physics to describe some aspects of social and political behavior 
From the theoretical point of view, opinion models are interesting to physicists because they 
present order-disorder transitions, scaling and universality, among other typical features of 


physical systems, which called the attention of many groups throughout the world 4l-ll5|. 


The basic ingredient of models of opinion dynamics is conformity, an important behavior 
of individuals that emerges as a result of their interactions with other individuals in the 
population [l, 2|. As examples: (i) an individual may copy the state (opinion) of one of 
his/her neighbors (the voter model 16l-ll9|). or (ii) he/she can consider the majority or the 


minority opinion inside a small group (the majority-rule models 20 


23|), or (hi) a given pair 


of individuals interact throught kinetic exchanges like an ideal gas 7|, [24,12^ . Among these 
models, we highlight the Galam’s majority-rule model 20|, l2l|. Indeed, influence of rnaiority 


opinions against minorities have been studied by social scientists since the 50’s 


30, 




However, recently the impact of nonconformity in opinion dynamics has attracted atten¬ 


tion of physicists 
independence [32 


26 


29| . There are two kinds of nonconformity, namely anticonformity and 


33| . and it is important to distinguish between them. The anticonformists 


are similar to conformists, since both take cognizance of the grou p n orm. Thus, conformists 


agree with the norm, anticonformists disagree. As discussed in 


32 


, l33| , an anticonformist 


actively rebels against influence. This is the case, for example, of the Galam’s contrarians 


2l| . individuals that known the opinion of the individuals in a group of discussion, and adopt 


the choice opposite to the prevailing choice of the others, whatever this choice is. On the 
other hand, we have the independent behavior. In this case, the agent also take cognizance 
of the group norm, but he/she decides to take one of the possible opinions independently 
of the majority or the minority opinion in the group ^ 31|- As stated by Willis in 321. 
“The completely independent person may happen to behave in ways which are prescribed or 
proscribed by the norms of his group, hut this is incidental. It should also be noted that pure 
anticonformity behavior, like pure conformity behavior, is pure dependent behavior”. 

In terms of the Statistical Physics of opinion dynamics, independence acts on an opinion 
model as a kind of stochastic driving that can lead the model to undergo a phase transition 
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28 Un fact, independence plays the role of a random noise similar to social temperature 
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Finally, another interesting and realistic kind of social behavior is usually called 


inflexibility. Individuals with such characteristic are averse to change their opinions, and 
the presence of those agents in the population affects considerably the opinion dynamics 


0 


34 


From the theoretical point of view, the introduction of inflexible agents works 


in the model as the introduction of a quenched disorder, due to the frozen character of the 
opinions of such agents. 


In this work we study the effects of conformity and nonconformity in opinion dynamics. 
For this purpose, we consider groups of 3 or 5 agents that can interact through the majority 
rule, but with the inclusion of disorder (inflexibility) and/or noise (independence). We 
analyze these effects separately in the standard majority-rule model, and all together, in 
order to study the critical behavior of the system induced by the mentioned effects. 

This work is organized as follows. In Section II we present separately in three subsections 
the microscopic rules that dehne the distinct formulations of the model, as well as the 
numerical results. These numerical results are connected with the analytical considerations 
presented in the Appendix. Finally, our conclusions are presented in section III. 


II. MODEL AND RESULTS 


Our model is based on the Galam’s majority-rule model 20l-l22|. We consider a fully- 


connected population of N = ha+ub agents with opinions Aor B concerning a given subject. 
In this sense, we are considering a mean-held-like limit, since each agent can interact with 
all others. In this case, the microscopic dynamics disregards correlations, that will be taken 
into account after, when we will consider the model on regular lattices. The opinions are 
represented by Ising-like variables Oj = ±1 {i = 1,2, and the initial concentration of 

each opinion is 0.5 (disordered state). We will consider three distinct mechanisms in the 
formulation of our model, namely majority-rule dynamics, inflexibility and independence. 
Our objective is to analyze the critical behavior of the system, and in this case we will 
consider separately in the following subsections three distinct cases: (i) the majority-rule 
model with independent behavior, (ii) the majority-rule model with inflexible agents, and 
(hi) the majority-rule model with inflexible and independent individuals. 
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A. Majority rule with Independence 


In this case, we consider that some individnals in the population can show a nonconformist 
behavior called independence 26l-l29|. The following microscopic rules govern the dynamics: 


1. A group of 3 agents, say (i, j, fc), is randomly chosen; 

2. With probability q all the three agents in the group will act independently of the 
opinions of the group’s individuals, i.e., independent of the majority/minority opinion 
inside the group. In this case, with probability / all the three agents flip their opinions 
and with probability 1 — / nothing occurs; 

3. On the other hand, with probability 1 — q the group follows the standard majority rule. 
In this case, all agents in the group follow the local majority opinion (if the opinion 
of one agent is different from the other two, the former flips alone). 


In the case where the 3 agents do not act independently, which occurs with probability 1 — 



model as a kind of stochastic driving that can lead the model to undergo a phase transition. 
In fact, independence plays the role of a random noise similar to social temperature jj]. 


We analyze the critical behavior of the system, in analogy to magnetic spin systems, by 
computing the order parameter 


O 






( 1 ) 


where (...) stands for time averages taken in the steady state. In addition to the time 
average, we have also considered conhgurational averages, i.e., averages over different re¬ 
alizations. The order parameter O is sensitive to the unbalance between the two distinct 
opinions, and it plays the role of the “magnetization per spin” in magnetic systems. In 
addition, we also consider the fluctuations y of the order parameter (or “susceptibility”) 


X = N{{0^)-{0Y) 


( 2 ) 
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FIG. 1. (Color online) Order parameter O versus the independence probability q for typical values 
of the flexibility /, for the mean-field formulation of the model (no lattice). The symbols correspond 
to numerical simulations for population size N = 10000 (averaged over 100 simulations) and the 
full lines represent the analytical prediction, Eq. (jl]). 


and the Binder cumulant t/, defined as [40 1 


U = 1 


3(02)2 • 


( 3 ) 


As we are considering a mean-field formulation of the model, one can follow Refs. 24l29| 
to derive analytically the behavior of the stationary order parameter. The behavior of O is 
given by (see Appendix 1) 

Mf 


0 = 1 - 


1-q 


1/2 


( 4 ) 


or in the usual form O ^ {q — qc)^, where 

Qc = qdf) = ’ ( 5 ) 

and we found a typical mean-held exponent /3 = 1/2, as expected due to the mean-held 
character of the model. The comparison of Eq. (jH) with the numerical simulations of the 
model is given in Fig. [H for typical values of the hexibility /. One can see an excellent 
agreement among the two results. Eq. ([5]) also predicts that there is an order-disorder 
transition for all values of / > 0, which was conhrmed numerically, see Fig. [2] (a). 

We also estimated the critical exponents for many values of /. As a typical example, we 
exhibit in Fig. [2] the hnite-size scaling (FSS) analysis for / = 0.7 (see pannels b, c and d). 
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FIG. 2. (Color online) Numerical results for the mean-field formulation of the model with inde¬ 
pendence. (a) Order parameter O versus q for typical values of / and population size N = 10000. 
It is also exhibited the hnite-size scaling analysis for / = 0.7 (pannels b, c and d). We obtained 
qc ~ 0.264, f5 K, 1/2, 7 « 1 and u Ri 2. Data are averaged over 100 simulations. 


The critical values qc were identified by the crossing of the Binder cumulant curves, as can 
be seen in the inset of Fig. |2] (b), and the critical exponents /?, 7 and u were found by the 
best collapse of data. For all values of / we found /5~l/2, 7 !^! and u ^ 2, which suggests 
a universality of the order-disorder phase transition. In particular, the numerical estimates 
of the exponent [3 agree with Eq. (|1]), that predicts /3 = 1/2 for all values of /. Notice that 
the exponents (3 and 7 are typical Ising mean-held exponents, which is not the case for u. 

I, and was 


This same discrepancy was observed in other discrete opinion models [2J, l29l. 


associated with a superior critical dimension Dc = 4, that leads to an effective exponent 
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FIG. 3. (Color online) Order parameter O versus the independence probability q for the 2D 
(triangular lattice, L = 100) and 3D (bcc lattice, L = 20) cases, considering / = 0.5 (a). One can 
see the typical behavior of a phase transition. We also shown O versus q for the model defined on 
a ID ring with N = 10000 sites (b). In this case, the results for distinct values of / suggest the 
absence of a phase transition. All data are averaged over 100 simulations. 


V = 1/2, obtained from u = Dci^' = 2. In this case, one can say that our model is in the 
same universality class of the kinetic exchange opinion models with two-agent interactions 


M, 


25 


29 


39| . as well as in the mean-held Ising universality class. 


To test the universality of the model under the presence of a topology, we simulated the 
dynamics on two distinct lattices, namely a two-dimensional triangular lattice and a three- 
dimensional body-centered cubic (bcc) lattice. In this case, the presence of a topology will 
introduce correlations in the system, and we expected that the mean-held results are not 
valid anymore. The lattices were built as folows. The triangular lattice was built from a hnite 
Lx L square lattice with extra bonds along one diagonal direction. In this case, each group 
of 3 agents is chosen as follows. First, we choose an agent at random, say i. Then, we choose 
at random two nearest neighbors of i (say j and k), in a way that each one of the 3 agents 
{i, j and k) is a neighbor of the other two agents, forming a triangle. On the other hand, the 
bcc lattice was built from a cubic structure with linear size 2L, and each group contains 5 
agents that were chosen as follows. First, we choose a random plaquette of 4 neighbor sites, 
forming a square. The hfth site is randomly chosen between the 2 possible sites in order to 
form a pyramid. A typical behavior of the order parameter as a function of q is shown in 















Fig. [3] (a) for both cases (2D and 3D, considering / = 0.5), where one can observe a typical 
behavior of an order-disorder transition. Considering distinct valnes of /, we performed 
a FSS analysis ir order to estimate the critical exponents (not shown). Thns, for the 2D 
lattice we obtained the same critical exponents of the 2D Ising model for all valnes of /, 
i.e., /9 ~ 0.125, v 1.0 and 7 ~ 1.75, and for the 3D lattice we obtained the same critical 
exponents of the 3D Ising model for all valnes of /, i.e., (3 ~ 0.32, u ^ 0.63 and 7 ~ 1.24 41 |. 
These resnlts snggest that considering a bidimensional (tridimensional) system the model 
is in the nniversality class of the 2D (3D) Ising model. Finally, we simnlated the model on 
an one-dimensional ring, where each 3-agents gronp was formed by a randomly chosen site 
and its two nearest neighbors. Typical resnlts for the order parameter as a fnnction of q 
are exhibited in Fig. |3] (b). In this case, the resnlts snggest that there is no order-disorder 
transition, as in the ID Ising model. In this case, considering the resnlts for ID, 2D and 
3D lattices and also for the mean-held case, one can say that the majority-rnle model with 
independent behavior is in the Ising model nniversality class. 

The comparative phase diagram (mean held x 2D x 3D) is exhibited in Fig. 01 where 
we plot for the fnlly-connected case the analytical solntion for qc, Eq. (jS]). The behavior is 
qnalitatively similar in all cases, snggesting a frontier of the form 


Qc = 


1 

1 + 0/ 


( 6 ) 


where a = 4 for the mean-held case, a ~ 33 for the triangnlar lattice and a ~ 6.5 for the 
bcc lattice (these last two valnes were obtained by a ht of the data exhibited in Fig. 0]). 

Notice that the mean-held analytical calcnlation, Eq. (jS]), overestimates the critical 
points qc, as it is common in mean-held approximations. However, onr calculations predict 
the occurrence of order-disorder phase transitions, as well as correctly predicts the form 
qc = 1/{1 + a f) between qc and /. 

From the phase diagram of this formulation of the model we can see that the increase of 
the hexibility parameter / leads to the decrease of qc, as also indicated in Eq. ([5]). This can 
be understood as follows. The increase of / leads the agents to perform more independent 
opinion changes or spin hips (which represents a nonconservative society). This action 
tends to disorder the system even for a small value of the independence probability q, which 
decrease the critical point qc. 
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FIG. 4. (Color online) Comparative phase diagram of the majority-rule model with independent 
behavior of the agents, for the mean-field, triangnlar lattice and bcc lattice cases, separating the 
Ordered and the Disordered phases. The symbols are the numerical estimates of the critical points, 
obtained from the crossing of the Binder cumulant cnrves for different population sizes. The lines 
are given by the analytical form, Eq. Q, with a = 4 for the mean-field case [full line, according 
to Eq. ([5])], a 33 (dashed line) for the triangnlar lattice and a ^ 6.5 (dotted line) for the bcc 
lattice, the last two values obtained from data hts. 


Notice that we obtained here for the mean-field case the same result for obtained in ref. 


26l |. where the independent behavior was considered in the Sznajd model. Indeed, in the 


mean-held formulation of the Sznajd model, the dynamics is very similar to the mean-held 
majority-rule dynamics for groups of size 3, which explains the identical result. However, in 
the mentioned reference, the model was not mapped in any universality class. 


B. Majority rule with Inflexibility 


As a second formulation of our model, we consider the majority-rule dynamics with the 
presence of some agents with the inflexibility characteristic, individuals whose stubbornness 


makes them reluctant to change their opinions [6|, 18|, 


34 


39] . As in [39|, we have considered 


a fraction d of agents that are averse to change their opinions. The following microscopic 
rules govern the dynamics: 
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FIG. 5. (Color online) Order parameter O versus the fraction d of inflexible individuals for the 
mean-field formulation of the model (main plot). The squares are the numerical results for pop¬ 
ulation size N = 10000 and the full line is the analytical prediction, Eq. ([7]). It is also exhibited 
in the inset the Binder cumulant curves for different sizes, showing the crossing of the curves for 
dc ~ 0.5, in agreement with Eqs. ([7]) and Q. Data are averaged over 100 simulations. 

1. A group of 3 agents, say {i,j, k), is randomly chosen; 

2. We verify if there is a majority of 2 (say i and j) in favor of a given opinion A or S, 
and in this case the other (say k) is a supporter of the minority opinion; 

3. If agent fc is a flexible individual, he/she will follow the local majority and flip his 
state Ok —)■ —Ok, otherwise nothing occurs. 


In this case, the frozen states of the inflexible agents work in the model as the introduction 
of a quenched disorder. As in magnetic systems one can expect that a disorder can 
induce/suppress a phase transition, as was also observed in the kinetic exchange opinion 
model with the presence of inflexibles 39j. 

As in the previous case (subsection A), one can derive analytically the behavior of the 
order parameter in this mean-held formulation of the model. The dependence of the order 
parameter with the the fraction d of inhexibles is given by (see Appendix 2) 


O = (1 -2d)^/2 ^ 


( 7 ) 
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FIG. 6. (Color online) Order parameter O versus the fraction d of inflexible individuals for the 
model with no independence defined on triangular lattices of distinct sizes L (main plot). It is also 
exhibited in the inset the Binder cumulant curves for different sizes, showing no crossing of the 
curves. Data are averaged over 100 simulations. 

or in the usual form O ~ (c? — dc)^, where 



and again we found a typical mean-held exponent jd = 1 / 2 , as expected due to the mean- 
held character of the model. The comparison of Eq. ([7]) with the numerical simulations of 
the model is given in Fig. [5l In addition, we also show in the inset of Fig. [5] the Binder 
cumulant curves for diherent population sizes, where one can observe the crossing of the 
curves at dc ~ 0.5, in agreement with Eqs. ([7]) and ([HD. Furthermore, a complete FSS 
analysis (not shown) give us/3«il/2, q^il and v k, 2, i.e., the same values obtained 
for the model presented in the subsection A. Thus, the presence of intransigents in the 
population leads the system to an order-disorder transition at a critical density d^ = 1 / 2 , 
and this transition is in the mean-held Ising model universality class. 

As a hnal observation of this subsection, we also simulated (as in the previous section) 
the majority-rule model with indexible agents on a two-dimensional triangular lattice, in 
order to test the universality of the model in comparison with the Ising model. The results 
are exhibited in Fig. [ 6 l One can see that the order parameter O (at least for the larger sizes) 
does not present the typical behavior of a phase transition, i.e., the usual change of concavity 
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of the curves. In addition, in the inset of Fig. |6]one can see that the Binder cumulant curves 
do not cross. A similar behavior was also reported in 39|] for another opinion model with 
inflexibility. In this case, those results suggest that the inclusion of inflexibility as we done 
here works in the model as a quenched disorder, and it destroys the phase transition in small 
dimensions like D = 2. In order to verify this hypothesis, we simulated the model on square 
lattices. In such case, we randomly choose a lattice site, and the group is formed by this 
random individual and his/her four nearest neighbors, forming a group of size 5 as in 43|. 
The behavior of O and U are very similar to the ones observed for the triangular lattice (not 
shown), suggesting that there is no order-disorder transition for D = 2. In some magnetic 


models such type of destruction due to quenched disorder was also observed [421. l44J. l45| 


C. Majority rule with Independence and Inflexibility 

As a third formulation of our model, we consider the majority-rule dynamics where agents 
can exhibit the independent behavior, as well as inflexibility. In this case, the model carries 
the rules of the two previous models (subsections A and B), namely: 

1. A group of 3 agents, say {i,j, k), is randomly chosen; 

2. With probability q each one of the three agents in the group will act independently of 
the opinions of the group’s individuals, provided he/she is not an inflexible individual. 
Thus, with probability / all flexible agents flip their opinions and with probability 
1 — / nothing occurs; 

3. On the other hand, with probability 1 — q the group follows the standard majority 
rule. In this case, each flexible agent follows the local majority opinion. 

Notice that, even if the agents decide to act independently of the group’s opinions, we 
will not see necessarily 3 changes of opinions, as in the model of subsection B. Indeed, the 
3 agents can change their opinions, but we can have two, one or even zero spin flips due to 
the frozen states of the inflexible agents. 

As in the previous cases, one can derive analytically the behavior of the order parameter 
as a function of the fraction d of inflexibles and the independence probability q, in the 
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FIG. 7. (Color online) Order parameter O versus the independence probability q for d = 0.2 
and typical values of the flexibility /, for the mean-held formulation of the model. The symbols 
correspond to numerical simulations for population size N = 10000 (averaged over 100 simulations) 
and the full lines represent the analytical prediction, Eq. ([9]). 

mean-field formulation of the model. The calculations give us (see Appendix 3) 

o = 1(1 - <i)2 (1 - (9) 


where 6 = 6{q, d, f) is given by 

e = 


''A 


( 10 ) 


( 11 ) 


Writing the order parameter in the the usual form O {q - Qcf, one obtains 

I-‘2d 

and again we found a typical mean-held exponent /9 = 1/2, as expected due to the mean-held 
character of the model. Notice that we recover the results of Eqs. ([5]) and (jS]) for d = 0 (no 
inhexibility) and Qc = 0 (no independence), respectively. The comparison of Eq. (j9]) with 
the numerical simulations of the model for d = 0.2 and typical values of / is given in Fig. 
[71 In addition, a complete ESS analysis (not shown) for many values of / give us (3 ~ 1/2, 
7 1 and u ^ 2, i.e., the same values obtained for the model presented in the previous 

subsections. Thus, this formulation of the model also leads the system to undergoes phase 
transitions in the same universality class of the mean-held Ising model. We also obtained 
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FIG. 8. (Color online) Phase diagram of the mean-field model with independence and inflexibility, 
in the plane q versus d for typical values of /. The symbols are the numerical estimates of the 
critical points, obtained from the crossing of the Binder cumulant curves for different population 
sizes. The full lines are the analytical form given by Eq. m- 

from the FSS analysis the critical points Qc for typical values of d and /. The comparison 
among the numerical estimates and the analytical prediction of Eq. mil is shown in the 
phase diagram of Fig. |8l 

From the phase diagram of this formulation of the model we can see that the decrease of 
the flexibility parameter /, related to the independent behavior, makes the ordered phase 
greater, for a given value of d < dc = 1/2. As in the case with no inflexibility, the increase of 
/ leads the agents to perform more independent spin flips, and this action tends to disorder 
the system even for a small value of the independence probability g, which decrease the 
critical point Qc- The presence of intransigent agents reinforces this behavior, leading to the 
decrease of the ordered phase for increasing values of d. 


III. FINAL REMARKS 

In this work, we have studied a discrete-state opinion model where each agent carries one 
of two possible opinions, ±1. For this purpose, we considered three distinct mechanisms to 
model the social behavior of the agents: majority-rule dynamics, inflexibility and indepen- 



15 


dence. Our target was to study the critical behavior of the opinion model under the presence 
of the mentioned mechanisms. Thus, we performed computer simulations of the model, and 
some analytical calculations complemented the numerical analysis. Let us remember that 
the original majority-rule model presents only absorbing consensus states with all opinions 


-|-1 or —1 


20 


22 |. 


First we considered the majority-rule model with independence in a fully-connected pop¬ 
ulation. In this case, there is a probability q that the 3 agents forming a group behave 
as independent individuals, changing opinion with probability / and keeping opinion with 
probability 1 — /. This mechanism acts in the system as a social temperature. In this case, 
we showed that independence induces an order-disorder transition for all values of / > 0, 
with the critical points Qc being a function of /. In addition, the model is in the same uni¬ 
versality class of the mean-held Ising model and of the kinetic exchange opinion models. In 
the ordered phase there is a coexistence of both opinions ±1, but one of them is a majority 
in the population. We observed that the larger the hexibility / (nonconservative societies), 
the smaller the value of q needed to disorder the system. In other words, for a population 
debating a subject with two distinct choices, it is easier to reach a hnal decision for a small 
hexibility concerning the independent behavior, as observed in conservative societies. Con¬ 
sensus states were obtained only for g = 0 or / = 0. As a test to the universality of the 
model, we simulated it on triangular and on bcc lattices, and we found the same critical 
exponents of the 2D and the 3D Ising models, respectively. In addition, simulations on ID 
lattices suggest the absence of a phase transition. All those results suggest the majority-rule 
model with independence is in the Ising model universality class. 

After that, we considered the majority-rule model in a fully-connected population with 
a fraction d of agents with the inhexibility characteristic. In this case, these agents present 
frozen states and cannot be persuaded to change opinion. In the language of magnetic 
systems, those special agents behave as the introduction of quenched disorder in the system. 
We showed that there is a critical fraction dc = 1/2 above which there is no order in the 
system, i.e., there is no decision or majority opinion. This order-disorder phase transition 
is also in the universality class of the mean-held Ising model. We observed consensus in 
the population only in the absence of indexible agents, i.e., for d = 0. In other words, the 
presence of intrasigents in the population makes the model more realistic, since there is a 
clear decision in the public debate for 0 < d < 1/2. Again, as a test for the universality of 
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the model, we simulated it on triangular and square lattices. In this case, we did not observe 
a phase transition for both lattices, suggesting that the model with quenched disorder does 
not undergo a phase transition in small dimensions like D = 2. 

Finally, we considered both effects, independence and inflexibility, in the majority-rule 
model. In this case we also observed a phase transition at mean-held level, and the critical 
points Qc depend on / and d. Consensus states were only obtained for g = d = 0, and the 
critical exponents are the same as observed before, i.e., we found again the universality class 
of the mean-held Ising model. From the phase diagram of this formulation of the model 
we observed that the increase of the hexibility parameter /, related to the independent 
behavior, makes the ordered phase smaller. 

It was recently discussed that the majority-rule model with limited persuasion can lead 


to the victory of the initial minority, provided it is sufficiently small 


46| . Thus, as a future 


extension of the present work, it may be interesting to analyze how diherent initial concen¬ 
trations of the opinions ahect the dynamics, as well as mechanisms of limited persuasion in 
the majority-rule dynamics. 


Appendix 


1. Analytical calculations: model with independence 


Let us consider the model with independent behavior in the mean-held formulation. Fol¬ 
lowing the approach of Ref. j^, 29|, we computed the stationary order parameter, as well 
as the critical values qdf)- Let us hrst dehne /i and /_i as the stationary probabilities 
of each possible state (-1-1 or —1, respectively). We have to calculate the probability that 
a given agent suffers the change -|-1 ^ —1 or —1 —)■ -|-1. We are considering groups of 3 
agents, so one can have distinct variations of the magnetization, depending on the states 
of the 3 agents. For example, the probability to choose at random 3 agents with opinions 
o = -|- 1 , i.e, a conhguration is ff. With probability 1 — q the conhguration re¬ 

mains (+,+,+), which does not affect the magnetization of the system. With probability 
g(l — /) the conhguration also remains and with probability qf the conhgura¬ 
tion changes to which cause a variation of —6 in the magnetization. In other 

words, the magnetization decreases 6 units with probability q f ff- One can denote this 
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probability as r(—6), i.e., the probability that the magnetization variation is equal to —6. 
Generalizing, one can dehne r{k), with —6 < k < +6 in this case, as the probability that 
the magnetization variation is k after the application of the models’ rules. As the order 
parameter (magnetization) stabilizes in the steady states, we have that its average variation 
must vanish in those steady states, namely, 

6 [r(+6) - r(-6)] + 2 [r(+2) - r(-2)] = 0 . (A.l) 


In this case, we have 


^^(+ 6 ) = fqf-i 
r{-G) = fqfi 

r(+2) = 3 (1 - g) /_i + 3 / g /i f\ 

r(-2) = 3 (1 - g) /i + 3 / g /_i . 

Thus, the null average variation condition Eq. flA.lj) give us 


(/i “ /-i) [/ q{fi + fi /-I + f-i) “ (1 “ q)fi /-I + f q fi f-i] — o (A. 2 ) 

which give us the solution fi = /_i (disordered state) or 

(1 - q) fi - (1 - 9 ) /i + /? = 0 , (A.3) 


where we used the normalization condition fi + /_i = 1. Eq. (]A.3|1 give us two solutions 
for /i, and the order parameter can be obtained from O = |/i — /-i|, which give us 


0 = 1 - 


4g/ 

1 - g 


1/2 


The critical points gc can be obtained by taking 0 = 0, 

1 

qc = 


1 + 4/ ’ 


(A.4) 


(A.5) 


that is the Eq. (jSD of the text. 


2. Analytical calculations: model with inflexibility 

Now we consider the model with inflexibility. Let us now denote the fraction of agents 
who have opinion +1 and are non-inflexibles by /i, and similarly for f_i. Notice that the 
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total fraction of inflexibles is d, and the fraction of inflexibles with opinion o = +1 is d/2, 
as well as the fraction of inflexibles with opinion o = — 1. In this case the normalization 
condiction becomes 

fi + /-I = I — d , (A.6) 

since the complementary fraction d represents the agents that have frozen states +1 or —1. 
Following the same approach of the previous appendix, the null average variation condition 
becomes 

2[r(+2)-r(-2)] = 0, (A.7) 

where the probabilities r{k) are given by 

r(+2) = 3 (/i + d/2)^ /_! 

r(-2) = 3 (/_! + d/2) Vi ■ 

Thus, the null average variation condition Eq. flA.Tp give us 


(/i-/-i)(/i/-i-c?V4) = 0, 

which give us the solution fi = f_i (disordered state) or 


(A.8) 


fi + {d — 1) fi + S '/4 — 0 , 


(A.9) 


where we used the normalization condition, Eq. flA.bp . Eq. flA.Qp give us two solutions for 
/i, and the order parameter can be obtained from O = \ fi — f-i\, which give us 


O = (1 -2d)^A 


(A.IO) 


The critical point dc can be obtained by taking 0 = 0, 

1 


4 - 2 > 


(A.ll) 


that is the Eq. (IHl) of the text. 


3. Analytical calculations: model with independence and inflexibility 


Now we consider the model with inflexibility and independence. As in the previous 
section, let us now denote the fraction of agents who have opinion +1 and are non-inflexibles 
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by /i, and similarly for /_i. Notice that the total fraction of inflexibles is d, and the fraction 
of inflexibles with opinion o = +lis(i/2, as well as the fraction of inflexibles with opinion 
o = — 1. The normalization condition is given by Eq. flA.611 

Following the same approach of the previons appendices, the nnll average variation con¬ 
dition becomes 

6 [r(+6) - r(-6)] + 4 [r(+4) - r(-4)] + 2 [r(+2) - r(-2)] = 0 , (A.12) 

where the probabilities r{k) are given by 

^^(+6) = fqf-i 
ri-G) = fqfi 

r{+4:) = Sfq {d/2) + 3 f q {d/2) f\ 

^^(-4) = 3fq {d/2) /i + 3 / g {d/2) f/ 

r{+2) = 3 f q (d/2)^ f_i -|- 3 (1 — g) (/i -|- d/2)^ /_i + 9 f q (d/2)^ f_i + 3 f q f\ fi 
r{-2) = 3fq {d/2)^ /, + 3 (1 - g) {f_, + d/2f /i + 9 / g {d/2f fi + 3fqf //_i 


Thns, the nnll average variation condition 


(/i - /-i) 1/9 (/? + /i /-I + /!i) + 2 /,<i(/, + /_i) + /+ /,/, /_, 

- (1 - 9)(/i/-i - ii^/4)l = 0 , 
which give ns the solntion fi = f_i (disordered state) or 


(A.13) 


{^ - q) fi + {^ - q) {d - 1) fi + [f q + {I - q) {d‘^/4)] = 0 , (A.14) 


where we nsed the normalization condition, Eq. flA.6D . Eq. flA.14D give ns two solntions for 
/i, and the order parameter can be obtained from O = \ fi — f-i\, which give ns 

0 = [(l-d)2(i_0)]i/2 ^ ^A.15) 


where 9 = 9{q, d, f) is given by 

9 = 


/g+ (1 - g) 


d^ 


(1 -g)(l-d)2 
The critical points qc can be obtained by taking 0 = 0, 

1 -2d 

qc = 


(A.16) 


(A. 17) 


l + 4/-2d 

that is the Eq. (fTTD of the text. Notice that we recover the resnlts flA.Sp and flA.lip for 
d = 0 (no inflexibility) and gc = 0 (no independence), respectively. 
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